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Conformational properties of regular dendrimers and more general hyperbranched polymer stars 
with Gaussian statistics for the spacer chains between branching points are revisited numerically. 
We investigate the scaling for asymptotically long chains especially for fractal dimensions df = 3 
(marginally compact) and d{ — 2.5 (diffusion limited aggregation). Power-law stars obtained by 
imposing the number of additional arms per generation are compared to truly self-similar stars. 
We discuss effects of weak excluded volume interactions and sketch the regime where the Gaussian 
approximation should hold in dense solutions and melts for sufficiently large spacer chains. 
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I. INTRODUCTION 


Hyperbranched stars with Gaussian chain statistics. 
Regular exponentially growing starburst dendrimers, as 
sketched in fig. ^ and more general starlike hyper¬ 
branched chains (l| with Gaussian chain statistics have 
been considered theoretically early in the literature [3-Q 
and have continued to attract attention up to the recent 
past I9I-II6I . One reason for this is that hyperbranched 
stars |17l.ll8l| with sufficiently large spacer chains between 
the branching points (as indicated by the filled circles) 
are expected to be of direct experimental relevance under 
melt or 0-solvent conditions |19l - [^ . Assuming a tree-like 
structure and translational invariance along the contour, 
the root-mean-square distance Rs between two monomers 
n and m, as shown in panel (a), is thus given by 

= {(Lm - Ln)^} = with ^ = 1/2 (1) 

being the inverse fractal dimension of the spacer chains, 
s the curvilinear distance along the tree between both 
monomers and b the statistical segment size of the spacer 
chains [l^. As a consequence, the typical distance Rg 
between the root monomer and the end monomers of the 
most outer generation g = G of spacer chains, as one 
possible observable measuring the star size [^, scales 
as = b'^SG with S being the length of the spacer 
chains (assumed to be monodisperse). Other moments 
are obtained from the normalized distribution P{r,s) of 
the distance r = |r„ — r„| which, irrespective of the 
specific topology of the branched structure, is given by 


P{r, s) 



( 2 ) 


with d = 3 being the spatial dimension [ 2 ^. Due to their 
theoretical simplicity such Gaussian chain stars (includ- 
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FIG. 1: Sketch of different topologies of branched polymer 
stars considered: (a) Regular dendrimer of generation num¬ 
ber G = 2 with M = 9 arms, (b) Hyperbranched so-called 
“a-star” with imposed spacer chain number Mg ~ g°‘~^ for 
g < G constructed iteratively (<?—>-( 7 - 1 - 1 ) by restricting ran¬ 
domly the branching of the arms. Some branches may thus 
end at (7 < G. (c) Self-similar fractal “/3-stars” are gener¬ 
ated starting with regular G = 2 dendrimers and replacing 
iteratively the Mi-i spacers of length Si-i by Mi = Mi-iUM 
spacers of length Si = Si-i/ns- The generator shown cor¬ 
responds to self-similar stars of constant density (df = 3). 
(d) Multifractal “ 7 -stars” are obtained by applying randomly 
more than one generator. Mixing with equal weight the gen¬ 
erator B {ns = um = 4) with the compact star generator A 
{ns = 4, um = 8) leads to a star with df = 2.5. 


ing systems with short-range interactions along the topo¬ 
logical network) allow to investigate several non-trivial 
conceptual and technical issu^ both for static [l3, Ell 
and dynamical @,[11EB [24 1^ properties, related 
to the in general intricate monomer connectivity imposed 
by the specific chemical reaction history. 

Aim of current study. We assume here that (i) the 
chemical reaction is irreversible (quenched), (ii) all 
spacer chains are monodisperse of length S and (Hi) flex¬ 
ible down to the monomer scale and (iv) that the branch- 
















2 


ing at the spacer ends is at most three-fold (/ = 3) as in 
the examples given in fig. [T] Our aim is to revisit vari¬ 
ous experimentally relevant conformational properties in 
the limit where the total monomer mass N and the to¬ 
tal number M = {N — 1)/S of spacer chains become 
sufficiently large to characterize the asymptotic univer¬ 
sal behavior and to sketch for different star architectures 
the regimes where the Gaussian spacer chain assump¬ 
tion becomes a reasonable approximation. We focus on 
the large-^ limit since this allows under 6>-solvent [2^ or 
melt conditions to broaden the experimentally meaning¬ 
ful range of the generation number G of spacer chains. 

Fractal dimension. One dimensionless property char¬ 
acterizing the star classes considered below is their fractal 
dimension d{ which may be defined as [l^, 


df = 


lim ‘"f") 

R^oo \og{R/b) 


( 3 ) 


with N being the mass and R the characteristic chain 
size. (Less formally, this definition is often written 
N ^ ( 111 -) For the regular dendrimers shown in 

panel (a) the number of spacers M and, hence, the to¬ 
tal mass N increase exponentially with the generation 
number G, while the typical chain size R{G) ~ VG only 
increases as a power law. That the fractal dimension 
thus must diverge, is denoted below by the shorthand 
“df = oo”. In addition we shall consider star classes 
of finite fractal dimension df, focusing especially on not 
too dense systems which should be (at least conception- 
ally) of experimental relevance. Specifically, we consider 
(i) marginally compact chains (29| of fractal dimension 
df = d = 3 and (ii) stars of fractal dimension df = 2.5 
which might be thought of as b eing assembled by diffu¬ 
sion limited aggregation (DLA) [24 . 

Power-law stars. As sketched in panel (b), such hy- 
perbranched stars of finite fractal dimension may be con¬ 
structed most readily by imposing a number of spacer 
chains Mg per generation g such that the power law 
Mg g°‘~^ holds. Hence, M N ^ G“. The “growth 
exponent” a of these so-called “a-stars” is set by the 
fractal dimension 


a = dfv (4) 

as may be seen using N ^ R'^^ and R ~ Re ^ {SGY 
[^ . While being a natural generalization of the regular 
dendrimer case, restricting the branching of star arms 
does, unfortunately, not lead to a self-similar tree since 
the iteration 5 —^ g -|- 1 is not a proper self-similar gen¬ 
erator acting on all spacer chains [^,1^. We therefore 
also consider truly self-similar (multi)fractal stars, called 
in the following /3- and 7-stars, generated iteratively as 
shown in panel (c) and panel (d) of fig. H] by the iterative 
application of a well-defined generator (or several gener¬ 
ators) on all the spacer chains as in the recent theoretical 
work on Vicsek fractals Q. For the latter architectures 
one thus expects to observe for the intramolecular coher¬ 
ent form factor F{q) the power-law scaling [l^ 1^ 

F{q) ^ for df < d = 3 (5) 


in the intermediate regime of the wavevector q. Note 
that eq. ([5]) only holds for open or marginally compact 
self-similar structures [2^|33|. In fact, Gaussian hyper- 
branched stars with higher fractal dimension, df > d, 
approach with increasing generation number and mass 
the Gaussian limit 

F{q) Ri A^exp {-{qRgf/d) for q < 1/bS^/^ (6) 

as shall be demonstrated below. 

Outline. The paper is organized as follows: We sum¬ 
marize first in sect. Ini the numerical methods and specify 
then in sect. m the different topologies studied. Some 
real space properties are presented in sect. IIVI before we 
turn to the characterization of the intramolecular form 
factor F{q) in sect. [V] While most of this study is ded¬ 
icated to strictly Gaussian hyperbranched stars, i.e. all 
excluded volume effects are switched off, we investigate 
more briefly in sect IVIl by means of Monte Garlo (MC) 
simulations effects of a weak excluded volume in¬ 
teraction penalizing too large densities. Even an expo¬ 
nentially small excluded volume is seen to change qual¬ 
itatively the behavior of large regular dendrimers. We 
conclude the paper in sect. I VIII Neglecting deliberately 
the long-range correlations expected as for linear chains 
[3^ . we sketch the regime where the Gaussian approxi¬ 
mation for melts of hyperbranched stars should remain 
reasonable for sufficiently large spacers. 


II. SOME COMPUTATIONAL DETAILS 

Settings and parameter choice. We suppose that the 
monomers are connected by ideal Gaussian springs. The 
spring constant is chosen such that the effective bond 
length b, eq. (HD, becomes unity. Also, both the tempera¬ 
ture T and Boltzmann’s constant are set to unity. All 
Gaussian spacers are of equal length S (which comprises 
one end monomer or branching monomer). With M be¬ 
ing the total number of spacer chains, a hyperbranched 
star thus consists of V = 1 -|- SM monomers. If nothing 
else is said. S' = 32 is assumed. (This arbitrary choice is 
motivated by simulations of dendrimer melts presented 
elsewhere.) For S = 32 we sampled up to a genera¬ 
tion number G — 11 for regular dendrimers and up to 
G r: 2000 for power-law hyperbranched stars of fractal 
dimension df = 3 and df = 2.5. (Even larger G obtained 
using smaller S are included below where appropriate.) 
Some properties of the largest system computed for each 
investigated star architecture are listed in Table [H 

Local and collective MC moves. Due to their Gaus¬ 
sian chain statistics many conformational properties can 
be readily obtained using Gaussian propagator tech- 
ni ques [l^ or equivalent linear algebra relations [ia& 
Il6l.l36j|. However, some interesting properties, such as the 
eigenvalues Af of the inertia tensor, can be more easily 
computed by direct simulation which are in any case nec¬ 
essary if long-range interactions between the monomers 
are switched on (see below). As shown in fig. [21 we use 
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star type 

d{ 

G 

A/10® 

Ae/10® 

(b 

r.f. 

Re 

Rg 

Dendrimer 

00 

17 

12.6 

197 

0.87 

0.11 

23 

22 

a-star 

6 

50 

22.6 

41.6 

0.74 

0.19 

40 

34 

a-star 

5 

80 

10.4 

10.0 

0.70 

0.22 

51 

42 

a-star 

4 

200 

7.2 

2.2 

0.63 

0.29 

80 

64 

a-star 

3 

2000 

16.2 

0.4 

0.47 

0.45 

253 

138 

a-star 

2.5 

2000 

2.4 

0.05 

0.36 

0.61 

253 

108 

(?-star 

3 

2048 

8.4 

1.2 

0.49 

0.51 

256 

179 

/3-star 

2.5 

4096 

1.1 

0.03 

0.45 

0.56 

362 

171 

7 -star 

2.5 

8192 

11.1 

0.3 

0.47 

0.54 

512 

351 


TABLE I: Various properties for different hyperbranched star 
types of spacer length S = 32: fractal dimension d{, largest 
generation number G, total mass N, number of end monomers 
Ne in the last generation shell g = G, rescaled Wiener in¬ 
dex (s)/smax with Smax = 2GS being the largest curvilin¬ 
ear distance between pairs of monomers, relative root mean- 
square fluctuation \/(s^} — (s)^/(s) (r-f-) of the normalized 
histogram w{s), root mean-square end distance Re between 
the root monomer and the end monomers of the generation 
shell g = G and radius of gyration i?g. 


(a) Pivot MC move (b) suitable data structure 



FIG. 2: Sketch of pivot MC move (a) and data structure 
(b). A monomer i (filled square) is selected randomly and all 
attached monomers k closer to the ends (within thin circles) 
are rigidly turned by an angle 0. A suitable data structure 
consists in ordering the spacer arms (their index indicated by 
the numbers) and the monomers such that all monomers k 
become neighbors in the monomer lists (i < k < j). 


pivot moves with rigid rotations of the dangling chain 
end (as shown by the monomers within the thin cir¬ 
cles) below a randomly chosen pivot monomer i. The 
monomers are collectively turned (using a quaternion ro¬ 
tation [1^) by a random angle 6 around an also randomly 
chosen rotation axis through the pivot monomer. As il¬ 
lustrated in panel (b) of fig. [H it is useful to organize the 
data structure such that arms and monomers which are 
turned together are also grouped together. This allows 
to rotate all monomer k with i < k < j. The tabulated 
monomer j = last(i), the last monomer to be turned, 
must be an end monomer. A pivot move does leave un¬ 
changed the distances between connected monomers. (If 
the connectivity of the monomers is the only interaction, 
a suggested move is thus always accepted.) To relax the 


local bond length distribution simple local MC jumps are 
added [33| . The root monomer at the origin never moves. 

Excluded volume interactions. Due to excluded vol¬ 
ume constraints the volume fraction occupied by a real¬ 
istic chain can, obviously, not exceed (much above) unity. 
One simple way to penalize too large densities is to in¬ 
troduce an excluded volume energy through the lattice 
Hamiltonian 

= l^n(r) (n(r) - 1) (7) 

r 

using the monomer occupation number n{r) of a simple 
cubic lattice. For all examples presented below we set 
Sx = 6y = Sz = 1, i.e. the grid volume 6V = Sx Sy Sz 
is unity and n(r) = p{r)6V measures the instantaneous 
local density. The Hamiltonian is similar to the finite ex¬ 
cluded volume bond-fluctuation model for polymer melts 
on the lattice described in [H, however, the particle 
positions are now off-lattice and only the interactions are 
described by the lattice. A local monomer or collective 
pivot move is accepted usin g th e standard Metropolis cri¬ 
terion for MC simulations [^. Note that the collective 
pivot moves are best implemented using a second lattice 
for the attempted moves. 

III. CHARACTERIZATION OF IMPOSED 
INTRACHAIN CONNECTIVITY 

Introduction. We assume that the hyperbranched 
star topology is not annealed, i.e. not in thermal equilib¬ 
rium, but irreversibly imposed by the chemical reaction. 
The first step for the understanding of such quenched 
structures is the specification and characterization of the 
assumed imposed connectivity, often referred to as “con¬ 
nectivity matrix” da, [ 3 . A central property charac¬ 
terizing the monomer connectivity is the normalized his¬ 
togram of curvilinear distances 

1 ^ 

( 8 ) 

n,m=l 

with Snm being the curvilinear distance between the 
monomers n and m. Trivially, w{s = 0) = 1/N and 
w{s) ~ 2N/N^ = 2/N for 0 < s ^ S' since the same 
monomer pair is counted twice. Note that the histogram 
w{s), sampled over all pairs of monomers of the chain, 
may differ in general from the similar distribution wo(s) 
of the curvilinear distances between the root monomer 
and other monomers. We remind also that for a linear 
polymer chain 

w(s) = —^ (1 -^ j for 0 < S < Smax (9) 

with Smax = — 1 ~ fV. For most of the star archi¬ 

tectures considered the largest curvilinear distance Smax 
is given by Smax = 2SG. The histogram w(s) will be 
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FIG. 3: Number of spacer chains Mg for dendrimers (bold 
solid line) and power-law stars of fractal dimension df = 3 
(a = 3/2). The open triangles refer to an a-star, the small 
filled triangles to a /3-star constructed as shown in fig. fTTc). 
The logarithmically averaged number of arms (thin line) and 
the root-mean-square fluctuations (circles) are of same order. 



FIG. 4: Number of monomers N/S ~ M vs. generation num¬ 
ber G for different imposed topologies and fractal dimensions 
d{. Dendrimers are indicated by df = oo (circles), a-stars by 
the other open symbols. The filled triangles corresponds to 
/3-stars of df = 3 {ns = 4, hm = 8 ) and df = 2.5 {ns = 16, 
u-M = 32) dimensions, the crosses to 7 -stars of df = 2.5. 

used below for the determination of experimentally rele¬ 
vant properties such as the radius of gyration Rg and the 
intramolecular form factor F{q). The first and second 
moments of ■u;(s) are given in Table |I] for the different ar¬ 
chitectures studied. We remind that N{s) is sometimes 
called “Wiener index” Wi S0. 

Regular dendrimers. Let us first summarize several 
simple properties of the regular dendrimers sketched in 
fig- HJa). As already mentioned above, the number Mg 
of spacer chains per generation shell g < G increases ex¬ 
ponentially as Mg = 3 • 2®“^ as shown by the bold line 
in fig. [31 Since we assume monodisperse spacer chains 
of length S, this implies wo{s) ~ 2'*/'® for S' ^ s < SG 



FIG. 5: Histogram w{s) measuring the number of monomer 
pairs at cnrvilinear distance s along the branched chain: (a) 
Dendrimers for G = 20 (bold line) and a-stars for different 
fractal dimensions (open symbols). The dashed line indicates 
the histogram for a linear chain of length N « Smax- Inset: 
Half-logarithmic representation for dendrimers. (b) Double 
logarithmic representation for df = 2.5 (upper data) and df = 
3 (lower data). As emphasized by the solid and dashed lines 
a power law x°‘~^ is observed only for j3- and 7 -stars. 

and that the mass N at generation number G must also 
increase exponentially, as shown in fig. j?) The histogram 
w{s) of curvilinear distances s for dendrimers is given in 
panel (a) of fig. |S] (bold solid lines). The main panel gives 
a linear representation of the dimensionless rescaled his¬ 
togram tc(s)smax as function of s/smax) the inset on the 
left-hand side a similar half-logarithmic representation. 
As one expects, the histogram increases exponentially 
for curvilinear distances S' <C s <C Smax due to the expo¬ 
nential increase of alternative paths of length s starting 
from an arbitrary monomer. Using simple combinatorics 
it can be seen that the histogram must become 

w{s) ^ ^ for 1 < S < Smax. (10) 

The cutoff observed for large s ~ Smax is due to the finite 
mass of the star and the finite length of its branches, just 
as the finite length of a linear chain gives rise to eq. ®. 
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As seen from Table HI the reduced first moment (s) /Smax 
approaches unity for dendrimers and the relative fluctu¬ 
ations are the smallest for all architectures considered. 

Hyperbranched a-stars. As already noted in the In¬ 
troduction, a simple way to generate stars of a finite frac¬ 
tal dimension d( is to impose a power law Mg = cg°‘~^ 
for the number of spacers per generation shell with c 
being a constant [s^. This is done by randomly attach¬ 
ing Mg spacer chains to the end monomers of generation 
g — 1 (with the constraint that at most two spacers can 
be attached per end monomer). An example for such an 
a-star with a = 3/2 is given in fig. [3] (open triangles). 
The corresponding total mass N Ri SG°‘ as a function 
of G is shown for a = 6/2 [s^, a = 5/2, a = 4/2 Q, 
a = 3/2 and a = 2.5/2 by open symbols in fig. S) The 
histogram wq{s) of curvilinear distances from the root 
monomer increases as wo(s) ^ for S' <C s < Smax 
as implied by the Mg-scaling (not shown). The curvilin¬ 
ear histograms w(s) over all pairs of monomers are pre¬ 
sented in the main panel of fig.[5Ka). The histograms are 
again non-monotonous increasing first due to the branch¬ 
ing and decreasing finally due to the finite length of the 
star arms. The latter decay becomes the more marked 
the weaker the branching, i.e. the smaller a, getting simi¬ 
lar for the smallest exponent a = 2.5/2 studied to the lin¬ 
ear chain behavior, eq. ®, indicated by the dashed line. 
As better seen from the double-logarithmic representa¬ 
tion in panel (b) of fig. [Sj a-stars cannot be describe d by 
a simple power law or exponential behavior for w{s) [40l| . 

Self-similar (3-stars. This is different for self-similar 
fractals created starting from a G = 2 dendrimer of 
spacer length Sq (as specified below) as initiator and iter¬ 
ating a generator as the one shown in fig. [TJc). At every 
iteration step i a spacer of length Si-i is replaced by 
um spacers of length Si = Si-i/ng. Hence, Si = So/n'g, 
Mi = 9n\f, Ni-1 = SiMi = 9So{nM/nsy and Gi = 2n^ 
for, respectively, the spacer length, the number of spac¬ 
ers, the total mass and the generation number of the 
star. Importantly, the arms added laterally to the origi¬ 
nal spacer can always be distributed such that the root- 
mean square end-to-end distance of the original spacer 
(filled circles) still characterizes the typical size of the re¬ 
placed spacer. Since SiGi = 2So for the curvilinear dis¬ 
tance between the root monomer and the end monomers 
in the largest generation shell g = Gi, the typical chain 
size R, thus remains by construction constant as we shall 
explicitly verify in sect. IIVI Note that the spacer length 
Si of the final iteration step is set by 

S = S, = So/n^s. (11) 

which fixes the mass A^o ~ <5'o ~ of the initiator star. 
Using Ni ^ ^ A'q* this implies 

um = ng with (3 = dfv (12) 

relating thus both numerical constants ns and um- As 
shown for df = 3 {ns = 2?, um = 2^) by the small filled 
triangles in fig. O such a self-similar construction leads to 


a strongly fluctuating number Mg of spacers. However, 
as shown by the thin solid line the (logarithmically) av¬ 
eraged number of arms still increases as Mg ^ g°‘~^ with 
a = /? = dfv in agreement with eq. ®. Interestingly, the 
corresponding (also logarithmically averaged) root-mean 
square fluctuations (as indicated by open circles) are of 
the same order, i.e. the relative fluctuations of spacer 
number Mg per generation shell are of order one. The 
important point is here that all monomers are statistically 
equivalent and that the root monomer does not play any 
specific role which would break the self-similarity. (As we 
have verified, this implies w{s) r; rco(s)-) Averaging over 
all spacer chains, the total mass N scales, as expected, 
again as N/S ^ G“ with a = dfV as shown in fig. IHby 
filled triangles for df = 3 and df = 2.5. The latter ar¬ 
chitecture, constructed using ns = 2'^ and um = 2 ^, is 
motivated by the fractal dimension df ~ 2.5 which may 
characterize self-similar stars generated by DLA in d = 3 
dimensions [1^ [30l - l^ . In our view this is one interest¬ 
ing universal limit of (at least conceptional) experimental 
relevance [2^ . Being self-similar all monomers are equiv¬ 
alent and since the number of monomers at a curvilinear 
distance s must increase on average as (s/S')““^, one ex¬ 
pects for S s Smax the power-law scaling 

^{S) ^ N X ^(s/5)“-l « ^(s/Smax)“-^ (13) 

^max 

with N r; 5'(smax/5')“ r; SG°‘. This is confirmed by the 
histograms (filled symbols) shown in fig. EJb). 

Stochastic two-generator multifractals. Since the 
DLA limit is of some importance we have sampled a 
second system class of fractal dimension df = 2.5 con¬ 
structed by mixing the generator A for marginally com¬ 
pact stars shown in panel (c) of fig. [T] with the second 
generator B shown in panel (d). Being constructed us¬ 
ing more than one generator these so-called “ 7 -stars” are 
in fact multifractals [ 2 ^, [s^ . (We remember that DLA 
clusters are also multifractal [^. No multifractal anal¬ 
ysis is required here, however.) For a given spacer 
we apply the generator A with a probability fA and the 
generator B with a probability fB = ^ — fA- By choos¬ 
ing different values of fA any fractal dimension between 
df = 2 and df = 3 can be sampled using both generators. 
By reworking the arguments leading to eq. (TT^l) it can be 
seen that /a = /b = 1/2 corresponds to df = 2.5. While 
/3-stars are deterministic, the 7 -stars have a stochastic 
topology due to the random mixing of both generators 
and an ensemble average over several stars is thus taken. 
As may be seen from the crosses in fig. [Hand fig.lHb), the 
properties of f3- and 7 -stars are, however, rather similar. 

IV. REAL SPACE CHARACTERIZATION 

End distance Rg. There are several ways to charac¬ 
terize the typical star size R all being equivalent from 
the scaling point of view. A double-logarithmic repre¬ 
sentation of the reduced mean-square end distance R^/S 
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FIG. 6 : Root-mean square end distance Re for different im¬ 
posed topologies: (a) Double-logarithmic representation of 
the reduced mean-squared end distance Re/S vs. reduced 
mass N/S. (b) Density p = N/R/. vs. total generation num¬ 
ber G for a spacer length S = 32. 


vs. the reduced mass N/S is presented in panel (a) of 
fig. ini Note that the values of Re obtained by direct 
MC simulations are within statistical accuracy identical 
to = b'^SG. Both data sets are lumped together. The 
regular dendrimer size increases, of course, logarithmi¬ 
cally with the mass (circles and bold solid line). The 
power-law slopes indicated for finite-df systems are con¬ 
sistent with the definition N ^ As one measure of 
the overall density of a star one may define p = N/Rf. 
(Obviously, a suitable order-one geometrical factor, such 
as 47 r/ 3 , might be included in this definition.) As can be 
seen from panel (b) of fig. [ 6 l the density for regular den- 
drimers exceeds already at G = 10 an unrealistic order of 
10 monomers per volume element. As indicated by the 
various power-law slopes, p ^ Q^-dv power-law stars 
of finite fractal dimension, i.e. the density increases for 
df > d and decreases for df < d as it should [ 1 , [ 2 ^ . 

Radius of gyration Rg. The radius of gyration Rg pre¬ 
sented in fig. [3 has been determined with identical results 
(lumped again together) either from the MC sampled 



FIG. 7: Reduced radius of gyration y = {Rg/Ref vs. genera¬ 
tion number G. The ratio y becomes constant only for /I-stars 
(filled symbols) and y-stars (crosses). 


configuration ensembles or by means the formula [M 


N 


r)2 _ _ 

s “ 2N^ 




n,m—l 


-| ^max 

- Y ( 14 ) 

s=0 


using the histogram of curvilinear distances w{s) dis¬ 
cussed above and the Gaussian chain property = b'^s. 
Measuring thus the first moment of w{s), the radius of 
gyration is equivalent for Gaussian chains to the Wiener 
index Wi. The reduced radius of gyration y = {Rg/Ref 
is plotted as a function of G. Since the end monomers 
dominate the mass distribution of regular dendrimers for 
large G, Rg becomes similar to Re- As expected, y ap¬ 
proaches unity from below (circles). Interestingly, the 
ratio y is constant for the self-similar /?- and 7 -stars, i.e. 
Re and Rg are similarly rescaled by the iterative applica¬ 
tion of the generators. This confirms the choice of gen¬ 
erators discussed in sect. uni We note finally that other 
observables characterizing R, such as the hydrodynamic 
radius [ 2 ^ , have been found to scale similarly as the end 
distance Re and the radius of gyration Rg. 

Density profiles. Figure[5]presents various normalized 
density profiles p{r) with r being the radial distance from 
the root monomer. The rescaled distribution y = p{r)R‘^ 
is plotted as a function of the reduced distance x = r/R 
with R = Re in panel (a) and R = Rg in panel (b) 
and panel (c). The distribution of the end monomers for 
regular dendrimers (G = 12, 5" = 32) shown in panel 
(a) is a reminder of eq. ([ 2 ]), i.e. of the trivial fact that 
the distances of all pairs of monomers have a Gaussian 
distribution (dashed line). The rescaled density p{r) = 
p{r)N of all monomers is shown in panel (b) of fig. [5] 
(using a half-logarithmic representation) for the largest 
star of each topology class. Note that the distribution 
p{r) has been either obtained for masses up to A Ri 10 ® 
from our MC simulations or for larger systems using 


Smax 

Pf") = Y '^o{s)P{r,s) 

s=0 


( 15 ) 
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FIG. 8: Density distributions y = p{r)R'^ with r/R being the 
reduced distance from the root monomer: (a) End monomer 
distribution with R — Re showing the expected Gaussianity 
(dashed line), (b) total monomer distribution rescaled with 
i? = TZg using the same symbols as in fig. [T] (c) double- 
logarithmic representation for three architectures with df = 
2.5. The slope indicates the exponent di — d = —0.5. 


with wo{s) being the already mentioned normalized his¬ 
togram of monomers of same curvilinear distance from 
the root monomer and P{r, s) the size distribution of a 
subchain of arc-length s given by eq. Since the den¬ 
sity distribution of large regular dendrimers (circles) is 
dominated by the end monomers, p{r) becomes essen¬ 
tially Gaussian (dashed line). We shall come back to this 
point at the end of sect. The histograms get natu¬ 
rally broader with decreasing df. Panel (c) on the right- 
hand side gives a double-logarithmic representation of 
the total monomer density distribution for three topolo- 

H i with df = 2.5. As explained in de Gennes’ book 
, the density should decrease as n{r)jr^ ^ ^j^d-df 
with n{j) ^ being the mass distributed within the 
volume r‘^. The same power-law exponent is obtained 
using u>o(s) ~ s““^ and integrating eq. (IT5)l for df < d 
and a; <C 1. Even the not self-similar a-star (open trian¬ 
gles) is seen to follow the predicted slope (solid lines). It 
is sufficient for this property that wo(s) has a power-law 
asymptotics albeit w(s) has not. 

Center of mass fluctuations. Albeit spherically aver¬ 
aged density profiles may reasonably characterize some 
aspects of the conformational properties of our hyper- 
branched polymer stars it is important to empha¬ 
size that a given instantaneous configuration may not 
be spherically symmetric and depending on the property 
probed experimentally or in a computer experiment these 
aspherical fluctuations become crucial. This issue is ad¬ 
dressed in fig. [HI The main panel compares the true ra¬ 
dius of gyration = T J2niiLn — Tcm)^) with a spher¬ 
ical approximation of the mass distribution defined by 
i?gp = assuming the center of mass of 

the star to be set by the root monomer at the origin 
for all configurations, i.e. = 0. The main panel of 



FIG. 9: Aspherical fluctuations: ratio {Rg/Repf (main panel) 
and rescaled largest eigenvalue A = (Ai) /R^ — 1/3 of the 
inertia tensor for regular dendrimers and a-stars (inset). 


fig- m presents (i?g/i?sp)^ as a function of G for different 
topologies. The ratio is always smaller than unity. The 
ratio is seen to approach unity from below for regular 
dendrimers and a-stars with df > d. While the spherical 
approximation = 0 becomes thus better with increas¬ 
ing size, stars with an incredible huge molecular mass 
are required to reach Rg m Rgp. Interestingly, the ratio 
decreases for a-stars with df = 3 and df = 2.5 (open tri¬ 
angles) while it is essentially constant for the self-similar 
(multi)fractals. For these experimentally most relevant 
star types the center-of-mass fluctuations remain thus 
relevant for asymptotically large chains. 

A sphericity. The asphericity of the stars may be 
(also) characterized by computing the three eigenvalues 
Ai > A 2 > A 3 of the inertia tensor of each star and aver¬ 
aging over the ensemble. Since R^ = (Ai) -I- (A2) + (A3), 
the rescaled eigenvalue A = (Ai)/i?g — 1/3 should vanish 
for perfectly spherical chains with (Ai) = (A2) = (A3). 
We have plotted A as a function of the inverse mass for 
several topologies in the inset of fig. [5] As expected from 
the consideration of i?sp, A is seen to vanish in the large- 
N limit for regular dendrimers and a-stars with df > d. 
(As shown by the solid line, A decays only logarithmi¬ 
cally with mass.) The opposite behavior is found for 
smaller fractal dimensions as shown by the open trian¬ 
gles. Whether for these systems A becomes constant for 
A —>■ 00 (as for linear chains) cannot be confirmed yet 
from our numerical data. 


V. FORM FACTOR 

Introduction. Gonformational properties of branched 
and hyperbranched star polymers can be determined ex¬ 
perimentally by means of light, small angle X-ray or neu¬ 
tron scattering experiments Using appropriate 

labeling techniques this allows to extract the coherent 


















intramolecular form factor F{q) defined as 


NF{q) = {p{q)p{-q)) 



exp(ig-r„) ||M (16) 


n—1 


with p{q) being the Fourier transform of the instanta¬ 
neous density and q the wavevector. The average is sam¬ 
pled over the ensemble of thermalized chains. For suffi¬ 
ciently large N and small q= ||q|| the radius of gyration 
Rg, as one measure of the star size, becomes the only 
relevant length scale. The form factor thus scales as [I^ 


F{q) = Nf{Q) with Q = qRg (17) 

being the reduced wavevector and f{Q) a universal scal¬ 
ing function with f{Q) = 1 — Q^/d in the “Guinier 
regime” for Q <C 1. The opposite large-g limit probes 
the density fluctuations within the spacer chains and the 
form factor becomes 


F{q) = 


12 


for 


6S'i/2 ^ ^ ^ 


1 

b' 


(18) 


For even larger wavevectors correlations on the monomer 
scale are probed. In the following we shall focus on the 
intermediate wavevector range l/Rg < q 1/bS^^^ be¬ 
tween the Guinier regime and the large-q limit. 

Dendrimers. Focusing on dendrimers, fig. llOl presents 
a Kratky representation of the form factor y = 
{F{q)/N)Q^ as a function of the reduced wavevector 
Q = qRg. Panel (a) shows stars of different spacer length 
S for a generation number G = 12, panel (b) different 
generation numbers G for a fixed spacer length S = 128. 
The increase of the rescaled data for very large wavevec¬ 
tors ( 7 ^ 1/6 observed in both panels is caused by the 
discrete monomeric units used in our simulations (see be¬ 
low) . The scaling observed for different S in panel (a) for 
the intermediate wavevector regime, where the Gaussian 
spacer chains are probed, is due to the fact that both the 
mass N and the radius of gyration R^ are linear in S. 
The corresponding failure of eq. (II3 in panel (b) shows 
that there is more than one characteristic length scale. 
Note that the strong decay after the Guinier regime above 
Q ~ 3 becomes systematically sharper with increasing 
generation number G. The bold solid lines in both pan¬ 
els indicate the expected asymptotic limit for G —> oo as 
discussed at the end of this section. Note that the den- 
drimer with G = 20 (large circles) shown in panel (b) is 
rather close to this limit. The form factor of this huge 
chain has not been obtained by MG simulations but by 
computing numerically the equivalent discrete sum 


^max 

F{q) = w{s)P{q,s) 

s=0 


(19) 


with w{s) being the curvilinear segment histogram dis¬ 
cussed above and P{q, s) the Fourier transform of the seg¬ 
ment size distribution P(r, s). Since for Gaussian chains 
P{q, s) = exp(—(aq)^s) with a = &/\/M, the form factor 




FIG. 10: Kratky representation of the form factor y = 
{F{q)/N)Q^ as a function of the reduced wavevector Q = qRg 
for dendrimers: (a) G — 12 for different spacer length S, (b) 
S = 128 for different generation number G. The dash-dotted 
lines indicate the Porod power law |^. the bold solid lines 
the predicted asymptotic behavior, eq. da¬ 


is readily computed yielding, as one expects, the same 
results as obtained from the explicitly computed config¬ 
uration ensembles. This can be seen from the dashed 
line in panel (a) of fig. [TU] for a spacer length S = 32. 
To compute numerically the form factor using w{s) has 
the advantage that the already mentioned discretization 
effect at q 1/b can be eliminated. To do this the dis¬ 
crete sum eq. (IT^ is replaced by a continuous integral 
for s > 0 and the s = 0 -contribution to the form factor is 
added. As shown by the thin solid line in panel (a), this 
allows to get rid of the irrelevant discretization effect. 

Marginally compact stars. Figure [TT] presents the 
form factor obtained using the continuous version of 
eq. (HU) for self-similar fractals of marginal compactness 
{df = 3). As one expects according to eq. ([5]), the data 
approach with increasing generation number the power- 
law slope 2 — df = — 1 (bold line) expected for the inter¬ 
mediate wavevector regime. We remind that eq. m can 
be derived from eq. (HU) and the scaling w{s) ^ ^ for 

self-similar fractals. Interestingly, eq. m does not hold 
for the (not self-similar) a-stars as may be seen from 
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FIG. 11: Kratky representation for /1-stars with df = 3. The 
reduced form factor approaches with increasing G the power- 
law slope —1 (bold line). The total monomer mass N is indi¬ 
cated for each iteration i. The dashed line has been obtained 
according to eq. (I20II by Fourier transformation of the spher¬ 
ically averaged density p{r) for i = 6. 



FIG. 12: Rescaled form factor y{Q) = {F{q)/N)Q^ for the 
largest stars available obtained using eq. dm). The self-similar 
P- and 7 -stars (filled symbols and crosses) decay, as expected, 
with a power law in the intermediate wavevector regime 

as shown by the thin solid line for df = 3 (i = 6) and by the 
dash-dotted line for df = 2.5. The dashed line indicates the 
preaverage approximation using eq. (Uni) for a-stars of df = 5, 
the bold solid line the expected large-dendrimer limit. 


fig- [l2l Note also that the large-g plateau of the rescaled 
form factor in fig. [TT] only decays as R^/N ^ 
extremely slowly with mass. This makes the numerical 
confirmation of the power-law slope demanding. For real 
experiments this implies that the determination of a frac¬ 
tal dimension df Ri 3 using the power-law scaling of the 
form factor for self-similar stars will also be challenging. 
We remind that a similar slow convergence of the inter¬ 
mediate wavevector regime is well-known for other more- 
or-less compact polymers such as polymers confined to 
ultrathin slits or melts of polymer rings [2^, . 


Comparison of different architectures. The rescaled 
form factors for the largest chains considered for each 
studied topology are compared in fig. [TH As expected, 
all data sets collapse in the Guinier regime below Q k, 2 
and become again constant for large wavevectors q 
1/65^/^. (The discretization effect for large q is again 
avoided using the continuous version of eq. (fTOl) .) The 
decay of the reduced form factor in the intermediate 
wavevector is seen to become systematically stronger 
with increasing fractal dimension df. For the self-similar 
stars this decay is described by eq. m as emphasized 
by the solid and the dash-dotted power-law slopes for, 
respectively, df = 3 and df = 2.5. All other architec¬ 
tures decay stronger than a power law. Note that it is 
the shape of this decay which is the most central prop¬ 
erty to be tested experimentally to characterize, at least 
approximatively, the structure of hyperbranched stars. 

Spherical preaveraging. As reminded at the beginning 
of this section, the intramolecular form factor is the en¬ 
semble average of the squared Fourier transform p{q) 
of the fluctuating instantaneous monomer density. Fol¬ 
lowing the recent work by Likos et al. [i^, this begs 
the question of whether in the limit of large and dense 
stars, where density fluctuations should become suffi¬ 
ciently small, one may replace p{q) by the Fourier trans¬ 
form p{q) of the averaged density profile p{r) discussed in 
sect. IIVl Due to the spherical symmetry of our stars this 
suggests using eq. (6.54) of ref. the approximation 

/ r sin((7 -r)\^ 

F{q) ij drp{r) ) (20) 

with p{r) = p{r)/N being known from eq. (1151) . As seen 
in fig. M eq. dini) is not useful for open (df < d) and 
marginally open (df ~ d) architectures for which the den¬ 
sity fluctuations are yet too large. The approximation 
becomes systematically more successful, however, with 
increasing fractal dimension as seen in fig. [T^ for a-stars 
of fractal dimension df = 5. Note that the striking decay 
of the rescaled form factor above the Guinier regime is 
accurately described by the approximation. As we have 
seen in hg. [H the distibution p{r) becomes systemati¬ 
cally more Gaussian with increasing star size and fractal 
dimension since the end monomers of the largest gener¬ 
ation shell dominate the total density. Since the Fourier 
transform of a Gaussian is again a Gaussian, this im¬ 
plies hnally eq. dH) as already stated in the Introduction. 
As seen by comparing the solid bold lines in fig. [TO] and 
fig. [TO] with the form factors computed using eq. (fTOl) for 
our largest dendrimers (circles), the asymptotic behavior 
eq. (I6|) gives an excellent fit to our numerical data. 

VI. WEAK EXCLUDED VOLUME EFFECTS 

Introduction. Up to now we have only considered ef¬ 
fects of the imposed monomer connectivity assuming all 
other interactions (persistence length, excluded volume, 
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FIG. 13: Radius of gyration Rg for dendrimers vs. excluded 
volume energy e for generation number G as indicated. Main 
panel: Unsealed raw data for S = 32. Inset: Data collapse 
of rescaled radius of gyration (Rg/i?*)^ as a function of the 
reduced excluded volume u(e)/u* with R* = Rg(e = 0) and 
Vi, = R'^/N^. The bold slope corresponds (approximately) to 
the compact limit N ~ Rg. 


...) to be switched off. Since essentially all proper¬ 
ties (apart the eigenvalues Xi of the inertia tensor) can 
be obtained analytically or numerically using the Gaus¬ 
sian chain statistics, the presented MC simulations were 
less crucial. Direct simulations are, however, essential for 
testing the influence of (albeit weak) excluded volume in¬ 
teractions computed using the lattice occupation number 
Hamiltonian, eq. described at the end of sect. m 
Scaling of chain sizes. Figure [T2] presents the ex¬ 
cluded volume dependence of the radius of gyration Rg 
for regular dendrimers. (Similar behavior is found for 
other characterizations of the typical chain size R.) As 
reveiled in the main panel, the excluded volume effects 
are the more marked the larger the mass N{G): The ra¬ 
dius of gyration increases already at e = 10“^ for G = 11 
while it has barely changed at e = 0.1 for G = 5. A suc¬ 
cessful data collapse is seen in the inset of fig. [T51 where 
the rescaled radius of gyration (i?g/i?*)^ is plotted as 
a function of the reduced excluded volume v{e)/vi, with 
Ri, = Rg{e = 0) ~ {SGY being the typical size of the 
Gaussian dendrimer star and 


;(e) = 5V (1 — exp(—/3e)) ~ jSeSV for /3e <C 1 (21) 


the excluded volume [20l| relevant for our model Hamilto¬ 
nian (/3 denoting the inverse temperature). The charac¬ 
teristic excluded volume u* below which the star should 
remain Gaussian is set by u* = R'f/N'^. This scaling is a 
direct consequence of Fixman’s general criterion 201 


1 » vp^Ri K. vN^/R^ 


( 22 ) 


for the Gaussian chain approximation with p Ri N/R'f the 
overall density for Gaussian stars. That the stars remain 
Gaussian for v/vi, <C 1 is emphasized by the horizontal 
asymptote indicated in the inset. The power-law slope 


7 = 1/d (bold line) for large reduced excluded volumes 
is only an approximative guide to the eye not taking into 
account logarithmic corrections. This can be seen (i) 
from the usual power-law ansatz [T^ Rg Ri Rf {v{e)/vi,y, 
(a) neglecting the weak logarithmic A^-dependence of i?* 
(fig. [6]) and (Hi) assuming that the dendrimers become 
essentially marginally compact, N ^ R'^, for large e in 
agreement with ref. [^. The latter point has explicitly 
been checked. For finite-df stars a similar scaling has 
been found (not shown). 

Spacer chain length criterion. We note finally that in 
terms of the generation number G and the spacer length 
5, Fixman’s criterion may be rewritten remembering that 
N Hi S 2^ for dendrimers and N ni S for power-law 
stars [1^ . Hence, the Gaussian approximation must hold 
for 5" <C S'* with an upper critical spacer length (2^ 

S* « {{v/bY 22G/Gdi-^-i/(2-dG 

S* « ((?;/6'^) (23) 

respectively, for dendrimers (df = oo) and finite-df hy- 
perbranched stars. In both cases S* Ri {b^/vY in 
d = 3 dimensions (while four-dimensional stars are only 
marginally swollen). 


VII. CONCLUSION 

Summary. We have revisited by means of direct an¬ 
alytical calculation, using for instance eq. (fTOl) . and MC 
simulations (sect.lU) several conformational properties of 
regular (exponentially growing) dendrimers and power- 
law hyperbranched stars (fig.© assuming Gaussian chain 
statistics = 1/2). As emphasized, a central imposed 
property is the normalized weight w{s) of curvilinear dis¬ 
tances s between monomer pairs (hg. [5]). Focusing on 
experimentally measurable observables such as the ra¬ 
dius of gyration Rg (fig. © and the intramolecular form 
factor F{q) (figs. [TUlfT^ . we investigated the scaling for 
asymptotically long stars with different fractal dimen¬ 
sions df. Due to their topological simplicity regular den¬ 
drimers (df = oo) have played a central role in our pre¬ 
sentation! fi£. [TOD as in other recent computational stud¬ 
ies [12, ka - fsil . Being (in our view) experimentally and 
technologically more relevant, we have also focused on 
stochastic architectures with df = 3 (marginally com¬ 
pact) and df = 2.5 as expected for stars created by DLA 
[^ . We compared “a-stars” constructed by imposing 
Mg ^ g°‘~^ arms per generation with truly self-similar 
so-called “/3-stars” and “y-stars” for which Mg becomes a 
strongly fluctuating quantity (fig.©- As shown in fig.dH 
only the latter two topologies show the power-law decay, 
eq. ®, of the form factor in the intermediate wavevector 
regime expected for open self-similar systems [12, IH, Ull . 
While large compact (df > d) stars may roughly be seen 
as dense colloidal spheres in agreement with Likos et al. 
[ 45 II , the instantaneous aspherical fluctuations cannot be 
neglected for experimentally relevant properties for the 
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FIG. 14: Sketch of critical spacer length S* for melts. The 
Gaussian star assumption holds above the bold lines. Note 
that the scaling argument does not allow to fix the scale of 
the vertical axis. If the generation number G is increased at 
constant spacer length S, as indicated by the dashed arrows, 
ideal chain behavior is expected for small G Gi,{S), while 
the star becomes colloid-like for larger G S> G*(S'). The 
number of chains interacting with a reference star should have 
a maximum at « G*(5'). 


smaller fractal dimensions studied (fig. Ill dashed line 
in fig. [HI)- We have commented briefly on the effects 
of gradually switching on an excluded volume potential. 
Coupling the (off-lattice) monomers by means of a (lat¬ 
tice) MC scheme (sect.lITl), we have sketched for different 
architectures the regime (e e*, S' <C S*) where the 
Gaussian star approximation can be assumed to be rea¬ 
sonable (fig. [HI). 

Conjectures for melts of hyperbranched stars. As al¬ 
ready pointed out, the Gaussian star assumption should 
be relevant under melt conditions assuming a large spacer 
length S ^ S*. That this holds can be seen by rewriting 
Fixman’s Gaussian chain criterion, eq. (1^^ . for melts 

1 » ^ p^Ri « N/Ri (24) 

remembering that the bare excluded volume v ^ e has 
to be rescaled by the total chain mass N [n, mm- 


The hyperbranched stars should thus remain Gaussian 
for interaction energies e ^ e* ~ k^TR'f /{N5V). Since 
e is not a parameter which can be readily tuned experi¬ 
mentally over several orders of magnitude, it is of some 
importance that eq. (1241) sets equivalently a lower bound 
Si, S depending on the generation number G. Follow¬ 
ing the discussion at the end of sect. I VII this implies 

Si, ~ {f2p for df = oo and 

Si, K inGt-d)u\ finite-df stars. (25) 

This scaling prediction is sketched in fig. |T4| for several 
architectures. Hyperbranched stars should remain thus 
Gaussian (albeit with a renormalized effective statisti¬ 
cal segment length [ 2 O, IH, js^) as long as ^ » S'*, if 
the interaction parameter fie is switched on as in the re¬ 
cent study of linear chain polymer melts [s^. Details 
may differ somewhat, of course, since the spacer chains 
may not be rigorously Gaussian due to long-range corre¬ 
lations related to the overall incompressibility of the melt 
. It is thus possible that even self-similar stars of im¬ 
posed df = 2.5 for the Gaussian reference (e = 0) may 
swell somewhat. We do conjecture, however, that this 
“swelling” for interacting large-S” hyperbranched stars in 
the melt remains perturbative as long as df < d = 3 [s^. 
Gonsidering the dynamical properties of strongly inter¬ 
penetrating hyperbranched stars for S Si, sampled us¬ 
ing standard molecular dynamics [s^, it will be of some 
interest to characterize the mean-square displacement of 
the star center of mass or, even better, the associated 
displacement correlation function [s^. As for the center 
of mass motion of linear polymer melts (ssl . , strong 
deviations from the Rouse scaling are to be expected [5^ . 
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